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Users  Guide  for  PNSEB 


a  program  for  rapid  computation  of  multioffset  vertical  seismic 
profile  synthetic  seismograms  for  layered  media 

INTRODUCTION 

This  document  describes  a  set  of  FORTRAN  programs  which  can  be  used  to  calculate 
synthetic  seismograms  and  a  few  other  response  characteristics  of  a  layered  elastic  medium. 
Copies  of  the  programs  can  be  obtained  fiom  D.  Lindwall  at  the  address  on  the  front  page 
of  diis  report  The  method  used  is  the  reflectivi^  method  as  described  in  Mallkk  and 
Blazer  (1987;  1988).  These  two  piq)ers  describe  the  background  of  the  metluxL  the 
calculations  used,  and  smne  of  the  li^tations  and  problems  invdved.  We  recommend  that 
the  user  cmisult  these  references  for  questimis  concerning  the  details  of  these  calculatXHis. 
Another  useful  reference  is  Rudman  et  aL  (1993)  whidi  reviews  the  two  Mallick  and 
Frazer  (1987;  1988)  papers,  gives  some  simple  «caiiq>les,  and  gives  the  complete 
FORTRAN  listing  of  the  basic  modules.  In  brief,  diis  set  of  programs  uses  a 
rearrangement  of  the  Kennett  reflectivity  algorithm  (Kennett  1983)  that  is  easily 
vectnizaUe.  The  Blon  method  of  quadrature  of  the  slowness  integral  is  used  for  an 
additKMial  increase  in  speed.  Coo^lex  fteqi^ncies  are  used  in  a  technique  to  avoid  time 
aliasing.  Attenuatkm  is  included  using  a  modified  Strick's  power  law  and  complex, 
finquency  dependent  velocities. 

An  earlier  version  of  this  code  was  for  modeling  OBS  (ocean  bottom  seismometer)  and 
OSS  (ocean  subbottom  seismometer)  data  and  was  called  SEABED.  This  eariy  code  was 
modified  so  that  it  could  model  the  very  long  range  Pq  and  Sq  waves  (or  Pq  and  Sq  in 
oceanic  crust)  and  was  renamed  to  PNSEB. 

Two  packages  are  used  sequentially  to  give  the  seismograms  or  other  desired  results. 
Alternative  versions  of  each  package  give  diffoient  types  of  results.  The  second  package 
produces  an  ASCII  file  containing  the  time  series  and  the  model  used  which  can  be 
transported  to  many  other  computers  and  used  for  making  plots.  Also  included  are  an 
ASCn  to  binary  conversicm  program  and  a  plotting  program. 

The  first  package  is  PNSN.  It  is  the  most  computaticmally  intoisive  package.  Using  a 
user  siq>plkd  modeL  it  calculates  the  earth  response  over  a  user  defined  finite  bandwidth 
for  a  ^tecified  range  of  receiver  offsets  and  receiver  depths  firom  a  single  source.  The 
calculation  is  done  in  fiequency-ray  parameter  (co-p)  space.  The  a>-p  response  is 
transformed  to  frequency-offset  (cn-x)  space  and  written  to  three  files  for  use  Ity  the 
seccxid  package.  The  pressure,  verticaL  and  horizontal  motitxis  are  each  written  to  a 
separate  file,  each  containing  the  response  for  all  receiver  depths  and  offsets.  A  version  of 
PNSN,  called  PNBES  uses  a  trapezoidal  integration  and  a  Bessel  function  instead  of  a 
Blon  integration  and  a  Hankel  function  as  in  PNSN  and  can  be  accurate  even  at  zero 
offset,  but  has  a  longer  run  time  than  PNSN  and  gives  poor  results  at  large  offsets.  Other 
current  versions  of  PNSN  calculate  only  the  pressure  response  (PNSN-P),  for  quicker 
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calculatiaii,  or  write  die  Cf>-p  response  (PNSN>WP)  so  that  it  can  be  transfraroed  to  tau-p 
by  a  versim  of  the  second  module. 

The  second  package  is  PNSYN.  It  reads  one  of  the  w-x  response  files  fimn  PNSN  or 
n^ES,  i^iplies  filters  and  windows,  convolves  the  response  with  a  source,  transforms  the 
response  to  the  time  dmnain  in  either  dme-horizontal  offset  (t-x)  or  time-depth  (t-z,  also 
calM  a  vertical  seismic  profile  or  VSP),  and  writes  an  ASCII  file  which  can  be  used  for 
plotting.  It  can  either  write  one  file  for  each  receiver  depth  or  one  file  for  each  receiver 
(^set,  depending  on  whether  t-x  seismograms  are  wanted  or  VSPs  are  wanted.  The 
version  SYNSPEC  produces  either  power  spectra  or  amplitude  spectra  of  the  response.  It 
includes  the  same  fflters,  windows  and  sources  as  PNS^.  The  version  TAUPSYN  uses 
the  output  from  PNSN-WP  to  generate  tau-p  synthetics. 

The  ASen  files  produced  by  PNSYN  and  SYNSPEC  can  be  transported  to  other 
computers  and  then  converted  CON  to  machine  dependent  binary  files  that  can  then  be 
used  by  VSPLT  to  produce  plot  files.  Some  users  have  developed  their  own  plotting 
routines  (which  may  use  commercial  graphics  software)  and  new  users  are  urged  to 
consult  with  old  users  to  learn  what  is  availaUe.  ROSECX)N  also  reads  the  ASCII  files 
from  PNSYN  and  SYNSPEC  and  writes  a  data  file  that  is  similar  to  ROSE  format 
(LaTraille  et  al.  1982;  LaTraille  and  Dorman  1983). 

PNSN  calculates  an  ipproximation  of  the  Green's  Function  (earth  response)  with 
numerous  limitations  caused  by  the  finite  sanpling  and  limited  bandwidth  in  both 
frequency  and  ray  parameter  space.  This  results  in  a  time  domain  "impulse"  response  of 
finite  width  that  always  has  sidelobes.  The  riupe  of  the  response  can  be  modified  by 
fiequency  and  ray  parameter  windows  and  frequency  filters.  Seismograms  of  the  response 
may  be  plotted  or  the  response  may  be  convolved  with  a  user-defined  source  wavelet 
Source  wavelets  are  usually  defined  to  be  cmvolved  with  a  delta  function  and  users  are 
cautioned  to  use  a  sufficiently  high  frequency  limit  in  the  response  calculation  and  a 
smooth  window  so  that  the  "deltaness"  of  the  calculated  response  is  sufficient  so  that  the 
source  wavelet  will  be  properly  reproduced.  An  additional  difficult  with  convolving 
source  wavelets  occurs  when  a  velocity  or  acceleration  response  is  calculated  rather  than  a 
displacemoit  response.  A  displaconent  "in^>ulse"  or  response  of  sufficiently  hi^ 
fi«quency  with  a  good  window  and  filters  will  approximate  a  delta  iiiq>ulse  response.  A 
velociQr  "in^ulse"  is  the  derivative  of  displacement  and,  in  order  to  follow  physical  laws, 
has  one  positive  pulse  and  one  negative  pulse  and  an  acceleration  "impulse,"  being  the 
derivative  of  velocity,  has  three  strong  pulses.  Measured  source  functions  for  velocity  and 
acceleration  devices  should  be  correct^  (first  or  sectmd  antiderivative)  so  that  they  can  be 
convolved  with  the  proper  response.  Geophones  goierally  measure  velocity,  or 
sometimes  acceleration,  and  the  proper  type  of  motion  must  be  calculated  so  that  the 
proper  amplitudes  are  obtained.  For  a  hydrophone,  even  though  it  measures  pressure,  or 
acceleration,  displacement  can  be  used  instead  to  get  a  simple  pulse  since  the  liquid 
medium  is  generally  uniform.  The  pressure  response  is  calibrate  to  give  amplitudes 
relative  to  the  response  at  a  distance  of  1  m  in  a  uniform  liquid  with  an  explosive  source 
(Ml  s  1,  M2  s  0,  M3  =  1  in  the  "resp"  file)  and  with  a  Hanning  window  (variable  IWIND 
s  2  in  the  "param"  file)  in  PNSYN.  The  velocity  and  acceleration  response  are  calibrated 
so  that  their  amplitudes  are  relative  to  the  response  at  1  m  below  a  vertical  inqiulse  source 
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fw  vertical  (M1,M2>I3»0,0,1)  and  1  m  hcffizontally  away  fire»n  a  horizontal  impulse 
source  for  horizontal  (Ml,M2>13sl,0,0)  with  a  Hanning  whidow  (vaiiaUe  IWIND  «  2). 
The  amplitudes  are  calibrated  for  the  Hannmg  window  since  it  is  the  most  useful  and  the 
least  fiequency  and  motion  type  dependent 

A  small  error  in  travel  times  occurs  when  using  very  low  Q  values  (high  attenuation) 
as  noted  Bromirsld  et  aL  (1992).  This  crrOT  comes  from  using  frequency  dependent  Q 
and  is  not  changed  by  different  values  of  the  WREF,  EPS,  and  SIG  variables  in  line  3  of 
the  "noodel"  file.  The  WREF.  EPS,  and  SIG  variables  do  affect  the  dispersion,  however 
and  the  user  may  wi^  to  select  your  own  values  rather  than  use  the  default  values, 
particulatly  when  and  are  very  low. 
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This  program  uses  the  reflectivity  technique  to  compute  the  syndtetic  seismic  response  for 
a  laterally  homogeneous  visco-elastic  ea^  and  may  be  used  to  model  both  e^loradon 
seismic  and  earthquake  data.  This  {xogram  is  eq)ecially  suitaUe  for  modeling  very  large 
offset  data  since  it  is  oqiaUe  of  handling  the  fiequencies  and  the  distances  invdved.  This 
program  is  written  for  use  on  a  Cray  X-MP  system  and  is  for  the  most  part  vectmized.  It 
has  been  run  on  a  Convex  con^uter  using  the  -cfc  option  in  the  FORTRAN  ocmqnler  and 
on  a  Sun  workstati<Hi  using  the  double  precisimi  q)tion  on  the  compiler.  We  have  found 
that  Sun  woricstations  using  the  IEEE  floating  point  standard  have  sufficient  accuracy  that 
they  can  be  run  with  single  precision  resulting  in  a  much  shorter  run  tinoe.  This  program  is 
also  ciqiaUe  of  computing  the  response  at  nuildple  receiver  depths  and,  thus,  may  be  used 
for  computing  multi-offset  synthetic  VSP  seismograms  as  well 

This  package  consists  of  the  main  program  and  six  subroudnes: 

1)  pnseb.f,  the  nudn  noodule. 

2)  The  subroutine  modlin.f,  used  to  read  the  model  file. 

3)  The  subroutine  respini,  used  to  read  the  file  "resp,"  which  craitains  computadonal 

parameters. 

4)  The  subroutine  executf,  used  to  compute  the  synthetics. 

5)  The  subroutine  bufoutf,  used  to  write  the  synthetics  to  a  disc  file. 

6)  The  logical  fiinctic.  cvmgtf. 

7)  The  subroutine  getsecs.f,  used  to  time  the  program. 

On  UNIX  systems,  the  con^iilation  and  linking  is  most  conveniently  done  with  a 
makefile.  A  sample  "makefile"  is  included  with  the  set  of  sample  input  files.  On  a  Ccmvex 
system,  the  -cfc  option  must  be  used  with  the  FORTRAN  compiler.  An  (ptimization  level 
of  2  (-02)  decreases  run  time  to  about  a  sixth  of  the  unoptimized  time.  Ccnpilation  time 
using  -02  is  about  8  minutes  on  a  Convex  Cl. 

This  program  can  use  a  large  anxnmt  of  menxnry.  On  some  machines,  the  charges  and 
the  swip  times  are  partially  dependent  on  the  memory  size  used.  To  minimize  the  memory 
required,  smne  of  the  array  dimensions  should  be  kept  to  the  minimum  required  for  each 
pr^leriL  In  the  subroutine  EXECUT,  the  memory  size  can  be  kept  down  by  pecifying 
array  dimensions  and  parameters  that  are  mily  as  Iv ge  as  necessary  for  the  calculation.  In 
the  main  commcm  block,  the  variables  CL,  QL,  CT,  QT,  RHO,  and  T  have  dimensions  of 
maximum  number  of  layers.  Be  sure  to  change  the  dimensions  for  these  vari^les  in  the 
COMMON  statements  in  all  of  the  other  subroutines  arxi  in  the  main  program.  The 
dimensions  for  the  variables  RMODEL  and  IRM  also  must  be  adjusted  for  the  maxiimim 
nuniber  of  layers.  The  parameters  NXMAX,  IBLK,  NZMAX,  IMAX,  and  UMAX  are 
used  as  the  dimensions  for  matty  cotiplex  arrays.  NXMAX  is  the  maximum  number  of 
(^Rtets  allowed.  It  can  be  set  to  one  for  calculating  a  single  VSP.  IBLK  is  a  blocking 
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parameter  to  reduce  the  memory  needed  by  the  program.  Program  speed  is  not  affected 
by  IBLK  unless  it  is  set  to  a  very  small  value.  N21MAX  is  set  to  one  when  there  is  only 
one  receiver  depth.  Using  multiple  receiver  depths,  as  in  the  case  of  a  VSP,  requires 
setting  NZMAX  to  a  number  equal  to  or  greater  than  the  number  of  receiver  depths  and 
will  turn  several  one  dimensional  con^lex  arrays  into  two  dimensional  coii^>lex  arrays. 
IMAX  is  set  equal  to  IBLK  *  NZMAX  and  ILMAX  is  set  equal  to  IBLK  *  the  number  of 
layers.  Every  time  that  any  of  the  array  dimensions  or  parameters  are  changed  in  the 
program  listing,  the  {Hogram  must  be  recompiled. 

The  maximum  nundwr  of  layers  the  program  can  haialle  is  adjustable  and  affects  the 
program  size.  This  is  changed  in  the  dimensions  of  the  arrays  ZZ,  CL,  CT,  QL,  QT, 
RHO,  T,  RMODEL,  and  IRM  in  modlinT  and  executf,  the  parameter  ILMAX  in  executf 
and  GL, . . .  ,T  in  the  transf.f  module  of  PNS  YN.  Most  of  these  arrays  have  one  element 
per  layer  and  those  that  don't  are  documented  in  the  source  listing.  Increasing  the  array 
dimensions  will  require  a  larger  memory  and,  on  some  systems,  will  slow  execution  time. 
Also,  be  sure  that  the  layer  array  dimensions  in  PNSYN  are  large  enough. 

INPUTS  TO  THE  PROGRAM 

The  inputs  to  the  program  are  in  two  files;  one  containing  the  model  and  one 
containing  the  program  parameters.  The  names  of  these  files  at  present  are  "model"  and 
"resp";  however,  these  names  can  be  changed  in  the  shell  script  "ronpn."  The  files  are 
written  in  free  format  with  either  commas  or  single  spaces  between  the  numbers.  In  these 
input  files,  the  names  of  integer  variables  begin  with  the  letters  I  through  N.  All  others 
variables  are  real  numbers  which  must  have  a  decimal  point 

The  "model"  file  contains  the  following  in  sequence: 

Line  1:  FS,  DECAY,  IFLT,  LSOUR,  IRND,  IQOPT  Where: 

FS  is  a  san^ling  parameter  for  the  ray  parameter,  usually  between  0.07  and  0.12. 
Sanq)ling  is  finer  as  FS  increases.  As  a  rule  of  thumb,  set  FS  at  0.07  when 
computing  synthetics  up  to  very  shcxt  ranges  Oess  than  100  km),  set  it  to  about 
0.09  for  computing  to  the  intermediate  ranges  Qess  than  1000  km)  and  set  it  to 
0.12  when  going  to  large  distances  (say  2000  km). 

DECAY  is  a  decay  factor  used  to  avoid  time  series  wraparound.  A  value  of  N  for 
DECAY  means  that  wrapped  around  signals  are  diminished  by  a  factor  of  N 
over  the  length  of  the  time  series.  A  rK:(xnmended  value  for  DECAY  is  SO.  If 
the  wraparound  is  in  the  negative  time  direction,  such  as  when  using  a  slow 
reduction  velocity,  use  a  DECAY  of  1  to  avoid  enlarging  the  wraparound. 

IFLT  is  the  earth-flattening  approximation  parameter.  If  IFLT=1,  then  the  program 
assumes  that  the  model  parameters  are  for  a  spherical  eartli  and  rq)plies  an 
earth-flattening  approximation  to  them.  (Xherwise,  the  model  parameters  are 
assumed  to  be  for  the  flat  earth  and  no  approximation  is  used. 

LSOUR  is  the  layer  on  top  of  which  the  source  is.  If  the  source  is  located  within 
a  layer,  then  divide  that  layer  into  two  layers  with  the  source  at  the  top  of  the 
second  layer. 


5 


IRND  is  a  flag  to  implement  nuidom  layer  velocities.  If  IRND  =  1.  then  line  2  is 
read. 

IQOPT  is  a  flag  to  reset  the  parameters  used  to  calculate  the  frequency  dq)eiKlent 
velocities  that  are  used  to  inoMrporate  realistic  attenuation  into  the  calculatkm. 
If  IQOPT  s  1,  then  line  3  is  read.  Otherwise,  the  default  values  are  used  for  all 
layers.  The  frequency-dependent  complex  velocities  are  calculated  with  a 
noodified  Stride's  power  law  (Mallick  and  Frazer  1987). 

Line  2:  (include  only  if  IRND  =  1;  otherwise,  delete  this  line)  LAYER,  SDLC,  SDCT, 

ZLAYER  where: 

LAYER  is  the  layer  number  for  the  tq>  layer  that  has  random  velodties.  All  layers 
(defined  in  line  4)  below  this  one,  except  for  the  bottom  layer  and  the  layer 
LSOUR,  have  velodties  that  are  randomized  over  a  unifonn  distribution  widiin 
the  specified  limits.  Random  values  are  added  to,  or  subtracted  from,  the 
values  given  for  the  parameters  CL  and  CT  in  the  line  4  list  The  parameters 
QL,  (^,  RHO,  and  T  are  unchanged. 

SDLC  is  the  standard  deviation  of  CL  for  the  random  velocities. 

SDCTT  is  the  standard  deviation  of  CT  for  the  random  velocities. 

2TAYER  is  the  thickness  of  the  random  velocity  layers  in  kilometers.  Any  layer 
thicker  than  this  is  divided  into  layers  that  are  ZLAYER  thick  or  less. 

Line  3:  (include  only  if  KJOPT  =  1;  otherwise,  delete  this  line)  LAYERN,  WREFP, 

WREFS,  EPSP,  I^SS,  SIGP,  SIGS.  These  are  values  for  the  variables  in  equation 

(28)  in  Mallick  and  Frazer  (1987)  and  are  used  to  calculate  the  frequency 

dependent  velocities  as  a  way  of  implimenting  realistic  attenuation. 

LAYERN  is  the  top  layo*  which  uses  the  user  supplied  values.  All  the  layers 
above  this  use  the  (tefault  values  and  all  the  layers  below  this  use  the  user 
tqrplied  values. 

WREFP  is  the  reference  frequency  for  compressional  waves.  The  default  value  is 
10  Hz. 

WREFS  is  the  reference  frequency  for  shear  waves.  The  default  value  is  10  Hz. 

EPSP  is  the  third  attenuation  parameter  for  compressional  waves.  The  default 
value  is  0.001. 

EPSS  is  the  titird  attenuation  parameter  for  shear  waves.  The  default  value  is 

0.001. 

SIGP  is  the  first  attenuation  parameter  for  compressicmal  waves.  The  default 
value  is  0.1. 

SIGS  is  the  first  attenuatitxi  parameter  for  shear  waves.  Thedefault  value  is  0.1. 

Line  4  (repeated):  CL,QL,CT,QT,RHO,T  (one  line  for  each  layer)  where: 

(X  is  the  compressional  wave  velocity  in  kilometers  per  second. 
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is  the  compiessuxial  Q. 

CT  is  die  diear  wave  velocity  in  IdlometoaAecond.  Note:  a  value  of  -1.  fm*  CT 
assumes  CT  «  CL/SQRT(3.). 

QT  is  die  shear  Q.  Note:  a  value  of -1.  for  QT  assumes  QT  »  QL/2. 

RHO  is  dw  drasity  in  grams  per  cubic  centimeter. 

Tis  the  absdute  diickness  cf  the  layer  in  kilometers. 

Note:  To  get  an  qifnoximadon  of  a  gradient  zone  between  two  layers,  the  user  may 
put  in  a  line  with  0 J^.T  for  die  first  three  variables  between  two  regular  layer  lines. 
This  will  automadcally  insert  N  layers  over  the  thickness  of  T  km  and  the  elastic 
values  will  be  linearly  interpolated  between  the  layer  above  and  below.  In  the 
gradient  ztme,  there  should  be  a  sufficient  number  of  layers  so  that  the  diickness  is 
equal  to  half  the  wavelengdi  of  the  highest  fiequency  (selected  in  "lesp"). 

The  re^XNise  file  "resp"  contains  the  following  in  sequence: 

Lines  1  and  2:  Utle  cards. 

Line  3:  NWJ^OR,TSECa*2W  where: 

NW  is  the  number  of  fiequencies  to  calculate.  The  maxinium  frequency  calculated 
will  be  NW  /  TSEC  hertz.  You  should  select  a  maxinium  frequency  about 
twice  that  of  what  you  want  to  (disove  since  a  harming  window  is  usually 
applied  over  fiequency  in  PNSYN  when  a  realistic  source  pulse  is  used.  If  a 
pulse  like  those  in  SAFARI  (Schmidt  1988)  is  used  (again  in  PNSYN)  be  sure 
to  calculate  to  a  fiequency  hi^  enough  for  that  pulse.  NWMAX  is  not  an 
array  dimensiai  parameter  in  executf  or  any  other  module  of  pnsn  but  it  is  in 
transf.f,  the  largest  module  of  PNSYN.  Be  sure  that  NTMAX  is  at  least  2  * 
NWMAX  in  transf.f. 

NOR  is  the  number  of  receiver  depths  for  a  vertical  seismic  profile  (VSP).  Set 
diis  to  1  to  get  a  horizontal  offset  (t-x)  seismograni.  If  more  than  one  depth  is 
used,  then  change  the  dimension  of  NZMAX  and  IMAX  in  executf  and 
NZMAX  in  transf.f  (one  of  the  PNSYN  modules)  so  that  they  have  the  same 
number  of  elements  in  the  rqiptopriate  dimension  as  there  are  receiver  depths. 
The  memory  space  needed  arid  output  file  sizes  are  neariy  proportional  to  the 
number  of  receiver  depths.  The  tun  time  is  only  slighdy  increased  v^ien  NOR 
is  increased. 

TSEC  is  die  time  series  length  in  seconds. 

P2W  is  the  maximum  ray  parameter  (p)  value  (in  seconds  per  kilometer)  over 
which  the  synthetics  ate  computed.  The  ray  parameter  is  the  inverse  of  the 
horiztmtal  compcment  of  the  {rfiase  velocity.  A  value  of  -1.0  for  P2W 
automatically  sets  P2W.  A  user-set  value  should  be  large  enough  so  that  the 
window  over  p  will  irx;lude  the  arrivals  with  the  slowest  horizontal  [diase 
velocity.  A  maximum  p  of  0.8  seems  to  be  sufficient  to  include  water  velocity 
arrivals. 
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line  4:  LOBS(I2S)»IZ>il»NDR  are  the  layers  on  top  n^icii  the  reodvers  are 
loaited.  If  receiver  is  witlm  a  layer,  introduce  a  fictidous  layer.  (The  source 
layer  is  qted:^  in  the  first  line  of  the  model  file.) 

Lines:  BX.FX,DX  where: 

BX  is  the  initial  range  in  kilometers. 

FX  is  the  final  range  in  Idlometers. 

DX  is  the  range  increment  in  Idkxneters. 

line  6:  PWl,  PW2.  PW3,  PW4 

These  are  ray  parameter  window  settings.  Responses  are  gently  tapered  between 
PW2  and  PWl  on  the  lower  side  (PW1<PW2)  and  between  PW3  and  PW4  on  the 
higher  side  (PW3<PW4)  using  cosine  banning  windows.  This  is  to  avoid  truncation 
louses  on  the  seisnoograms.  If  PWl.  PW2.  PW3,  and  PW4  are  all  set  to  zero,  then 
the  default  ti^  widths  of  15%  are  used  the  program  which  work  quite  well  If 
PW2  is  set  equal  to  PWl  or  PW3  equal  to  PW4.  then  the  default  values  are  used 
to  avoid  numerical  errors.  If  the  user  wants  to  see  reflections  whidi  have  a  anall 
ray  parameter  (such  as  near  vertical  iiKidence),  dien  you  may  need  to  use  a  small 
value  for  P\ll^  such  as  0.1, 0.05  or  even  0.01. 

line?:  DPMAX 

The  largest  step  size  in  p  used.  The  p  step  size  is  set  automatically  as  a  function  of 
frequency  except  when  DPMAX  is  smaller  than  the  automatic  value,  then  DPMAX 
is  used.  This  avoids  very  large  steps  in  p  at  low  frequencies  and  is  used  when 
luniung  PNBES  since  the  automatic  value  is  much  too  big  in  a  tr^zmdal 
integration. 

line  8:  ETYPE 

Set  this  source  parameter  to  D  if  the  source  input  are  straight  moment  tensor  and 
to  E  if  the  moment  tensra*  is  to  be  computed  the  earthquake  focal  mechanism 
solutkms.  PCX'  an  explosive  source  in  water,  use  D. 

Line  9:  (H(I),  1=1,2) 

This  is  the  linear  part  of  the  source.  H(l)  denotes  the  vertical  and  H(2)  denotes 
the  horizraital  part  For  earthquake  sources,  both  H(l)  and  H(2)  are  usually  zero. 
Fot  an  explosive  source  in  water,  set  H(l)  and  H(2)  to  zero. 

Line  10:  MO,  DELTA,  LAMDA,  PHIS,  PHI  or  Ml,  M2,  M3 

This  can  be  of  two  types.  If  line  8  is  E  then  line  10  contains  MO,  DELTA, 
LAMDA,  PHIS,  PHI,  where  MO  is  the  seismic  moment  DELTA  is  the  dip, 
LAMDA  is  the  rake,  and  PHIS  is  the  azimuth  of  the  fault  plane  and  PHI  is  tte 
azimuth  of  the  receiver  location.  If  line  8  is  D,  then  line  10  contains  Ml,  M2,  M3, 
where  Ml,  M2,  M3  are  the  (1,1),  (1,2),  and  (2,2)  components  of  the  moment 
tensor.  This  is  a  four-conqxment  2-D  renscx  instead  of  the  nine-component  3-D 
tenscx^that  you  will  find  in  textbooks.  The  (2,1)  component  is  the  same  as  the 
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(1^)  conqxxient  For  an  explosive  source  in  water,  set  Ml  and  M?  to  1.0  and  set 
M2  to  zero  (0.0). 

OITTPUT  OF  THE  PROGRAM 

The  program  has  four  output  files.  The  first  file  is  named  "pnsnout,"  and  contains  the 
conq}utati<Mi  information,  e.g.,  model,  ray  parameters  used,  total  CPU  time  used  m 
computing  etc.  This  is  a  text  file  and  may  be  written  to  the  screen  to  verify  whether  the 
job  has  executed  successfully.  The  other  files  are  named  “denp,”  "denh,"  and  "denv"  and 
omtain  the  computed  pressure  and  the  horizontal  and  the  vertical  responses,  req)ectively, 
in  frequency  range  domain.  The  amplitude  of  the  response  files  are  relative  to  the 
respcmse  at  1  m.  These  files  are  used  as  input  to  the  package  PNSYN  for  generating  the 
synthetics.  The  names  of  these  output  files  can  be  chwged  by  modifying  the  three  OPEN 
statements  iq>pearing  near  the  beginning  of  executf. 

RUNNING  THE  VERSION  PNBES 

PNBES  uses  exacdy  the  same  format  for  the  "model"  and  "respin"  files.  In  fact,  all  of 
the  subroutines  are  the  same  except  for  executf,  which  becomes  exebes.f,  and  a  new 
subroutine  called  bessub.f,  which  contains  the  Bessel  function.  A  few  of  the  parameters  in 
the  "resp"  file  must  be  changed  both  to  take  advantage  of  and  to  consider  sane 
restrictions  in  the  use  of  the  Bessel  function.  Itie  trapezoidal  ray  parameter  integration  in 
PNBES,  which  uses  a  polynomial  expansion  of  a  B^sel  function,  is  valid  even  when  the 
product  of  the  fiequency  (w),  ray  parameter  (p),  and  offset  (x)  (w*p*x)  is  small  while  the 
Filon  integration  in  PNSEB,  which  uses  an  asymptotic  expansion  of  a  Hankel  function  is 
invalid  for  a  small  w*p*x  than  with  the  F^on  integration  in  PNBES.  The  big  drawback 
with  the  trapezoidal  integration  in  PNBES  is  that  many  more  integration  points  are 
needed,  especially  for  large  w*p*x.  This  makes  PNBES  useful  for  calculations  when 
w*p*x  is  small  but  impractical  when  w^p'^'x  is  large. 

The  window  over  ray  parameter  (p)  must  be  changed  fiom  the  default  values  to  get 
the  correct  amplitudes  from  the  arrivals  with  a  small  w*p*'x.  The  default  p  window  is  a 
15%  cosine  taper  on  both  sides.  This  minimizes  the  artifacts  but  also  eliminates  any  real 
arrivals  with  a  small  p.  The  variables  PWl,  PW2,  PW3,  and  PW4  on  line  6  of  the  "resp" 
file  control  the  p  window.  Values  of  0.0, 0.01, 0.7,  and  0.8  with  an  oceanic  crustal  model 
give  a  good  comprcnnise  between  minimizing  the  truncadon  phase  and  keeping  most  of 
die  amplitude  at  small  x.  A  value  of  0.0  cannot  be  used  for  PW2  because  it  result  in  a 
numerical  error  and  halt  the  program.  A  very  small  PW2  will  give  truncadon  phases  with 
an  infinite  phase  velocity  from  the  stronger  arrivals  appearing  at  x  =  0  Im.  These 
truncation  phases  are  very  small  when  PW2  is  0.1  or  larger. 

The  parameter  DPMAX  on  line  7  of  the  "resp"  file  must  be  sufficiendy  small  so  that 
the  trapezoidal  integration  has  sufficient  sampling.  A  value  of  0.001  or  smaller  may  be 
needed.  Decreasing  DPMAX  increases  the  run  time.  TSEC  (line  3)  should  be  at  least 
twice  as  loig  as  the  desired  time  section  to  avoid  some  strong,  time-aliased  artifacts, 
which  appear  at  time  TSEC-T  in  the  record  where  T  is  the  true  arrival  time.  There  are 
also  some  very  strong  artifacts  in  the  last  2  s  of  the  time  series.  Remember  to  increase 
NW  proportional  to  TSEC  so  that  the  maximum  frequency  calculated  will  remain 
constant  Increasing  NW  also  increases  the  run  time. 
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SAMPLE  FILES: 

Following  tfiis  paragraph  are.  in  <»der.  a  "makefile"  fw  PNSN,  a  shell  script  file  to  run 
PNSN,  a  "model"  file,  and  a  "resp"  file.  For  PNBES,  change  execut.o  to  exebes.o  and 
change  pnsn.02  to  pnbes.Q2  in  the  makefile  and  change  pnsn.02  to  pnbes.Q2  in  the  shell 
script  file.  PNBES  uses  the  same  "model”  and  "resp"  files  as  PNSN. 


#  Makefile  for  pnsn 

#  Makefiles  are  used  In  UNIX  based  systems  and  are 

#  system  dependent.  The  best  way  for  you  to  make  such  a 

#  file  Is  with  the  "makemake”  command  and  then  edit  It  to 

#  Include  the  options  that  you  want.  This  makefile  Is  for 

#  a  CONVEX  system  and  includes  some  vectorization  and 

#  optimization  options. 

OBJ  =\ 

pnseb.o  modlin.o  respin. o  execut.o  bufout.o  cvmgt.o 
getsec.o 

FC  =  fc  -vn 

OPT  =  -01  -cfc 

PROF  = 

LINKPROF  = 

LIB  =  -Iveclib 

.  f  .o: 

$(FC)  $(0PT)  $(PR0F)  -c  $A.f 

pnsn.02:  $  (OBJ) 

$(FC)  $  (LINKPROF)  -02  -o  pnsn.02  $ (OBJ)  $  (LIB) 

#  cvmgt.fil:  cvmgt.f 

#  $(FC)  $(OPT)  -11  cvmgt.f 

execut . o : 

$ (FC)  $ (PROF)  -02  -nip  -cfc  -c  execut. f 
getsecs.o: 

$(FC)  $(PROF)  -c  getsecs.f 
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The  next  file  is  a  shell  script  (UNIX)  for  running  the  piognunPNSN.  It  assigns  logical 
numbers  to  files  and  directs  error  listing  to  "pnsn.ouL" 


#  model  is  the  model  Input  file 

#  resp  Is  the  response  parameter  file 
setenv  FOROlO  /lo2way/lindwall/pn/pnsn/model 
setenv  FOROll  /io2way/lindwall/pn/pnsn/resp 
setenv  FOR013  /io2way/lindwall/pn/pnsnout 
/lo2way/llndwall/pn/pnsn/pnsn.02  >&•  pnsn.out 

#  the  output  files  with  the  w-p  domain  synthetics  are: 

#  denp  with  the  pressure 

#  denv  with  the  vertical 

#  denh  with  the  horizontal 

The  next  file  is  a  sample  model  file  for  PNSN.  The  first  line  contains  some 
parameters.  The  second  line  is  the  first  model  layer,  10  m  diick.  The  receivers  are  at  the 
t(^  of  the  sectMid  layer  (variable  lobs  set  in  the  "resp”  file  below)  and  the  source  is  at  dte 
tq)  of  the  third  layer  as  set  by  the  variable  LSOUR  in  the  first  line  of  this  file.  Note  that 
the  first  three  layers  have  identical  physical  pr(q>erties  and  are  specified  as  different  layers, 
since  the  receivers  and  source  are  at  the  of  layers  two  and  three,  respectively.  The 
seventh  line  ^lecifies  a  "gradient,"  vriiidi  is  i^^proximaied  by  12  diin  layers  with  a  total 
diickness  of  1.8  km.  The  last  laya  is  set  to  be  duck  to  avoid  reflecticmi  finun  its  bottom 
but  a  default  built  into  the  code  makes  it  much  thiclar  than  this. 

0.07, 50. ,0,3, 0,0 

1.5.5000. .  0.. 0. . 1.03. 0.01 

1. 5. 5000..  0..0.. 1.03. 0.02 

1.5.5000.. 0..0..1.03.4.47 
2. 0,100. ,0.7, 10. ,1.5, 0.22 

5. 3. 450..  2. 8. 100.. 2. 4. 0.7 

0., 12., 1.8 

6.8. 1000..  3.7. 600.. 3.3. 0.6 

7. 5. 1000..  4. 05. 1000.. 3. 0.0.1 

8. 3. 1000..  4. 56. 1000.. 3. 3. 20.0 
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The  next  file  is  a  sanqile  "leq)"  file  for  running  tte  program  PNSN.  The  10-s  long 
time  series  and  512  fiequc^es  mean  that  tlw  inogram  will  calculate  the  response  at  every 
0.0195  Hz  (10/512)  up  to  a  maximum  of  51.2  Hz. 


; MODEL  CRUST 1 
:PVH 

512,1,10.0,0.8 

2 

5. 0,20. ,5.0 
0., 0.1, 0.7, 0.8 
0.002 
D 

0.,0.,0. 

1., 0.0,1. 


!nw, nor, t sec, p2w 
ilobs (nor) 

!bx,  fx,clx 

!  pwl , pw2 , pw3 , pw4 

!  dpmax 

!etype 

!h(l)  ,h(2) 

!ml,m2,m3 
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PNSYN 


This  program  is  to  be  used  to  process  output  from  the  package  PNSN/PNBES.  PNSYN 
reads  one  die  frequency-ray  parameter  domain  (w-p)  output  files  (such  as  denp)  of 
PNSEB,  filters  and  windows  dw  data,  involves  it  with  a  source  function,  adds  nmse,  and 
transformes  it  to  the  time  domain,  generating  the  synthetic  seismograms.  This  program  is 
written  for  running  on  a  Gray  X-MP  system  and  modified  to  run  on  a  CONVEX  using  the 
-cfcqition.  The  program  consists  of  the  main  program  and  eight  subroutines: 

1)  pnsyn.f,  the  main  program 

2)  subroutine  pannsin.f,  reads  the  parameter  file. 

3)  subroutine  trBnsf.f,  generates  synthetics. 

4)  sutvoutine  filter.f,  filters  the  data. 

5)  subroutine  fftf,  does  a  fast  Fourier  tiansfrmn. 

6)  subroutine  bufin.f,  reads  die  unformatted  (binary)  files. 

7)  subroutine  sourceS.f,  bubble  source  frinction. 

8)  subroutine  pulseLf,  SAFARI  pulse  type  1. 

9)  subroutine  getsecs.f,  times  the  program. 

10)  subroutine  makenoise.f,  makes  noise  then  shtqies  it  by  filtering. 

11)  subroutine  filtem.f,  filters  the  noise. 

12)  sulxoutine  ran2.f,  a  random  number  ^nerator  from  "Numerical  Recepies." 

13)  subroutine  gasdev.f,  turns  "ran2"  into  a  Gaussian  (normal)  distribution. 

INPUTS  TO  THE  PROGRAM 

There  are  two  input  files  to  this  program: 

1)  The  parameter  file  "param."  This  file  can  be  written  as  a  free  format  file,  whoe  the 

variable  type  is  the  default  taken  from  the  variable  name  with  integer  names 
starting  with  the  letters  I  through  N  and  the  rest  are  real  numbers.  Real  numbers 
must  have  a  decimal  point  but  integers  must  not.  The  name  of  diis  file  is  assumed 
to  be  "param”  but  the  name  can  be  changed  in  die  shell  script  "ninsyn.” 

2)  The  response  file,  which  is  the  output  from  PNSN,  is  one  of  three  files.  The 
pressure  respcnisels  "denp,"  the  horizontal  displacement  is  "denh,"  and  the  vertical 
displacemoit  is  "denv."  The  file  to  be  transformed  must  be  renamed  to  "synin"  (or 
copied  to  "synin").  The  name  of  the  input  file  is  assumed  to  be  "synin”  but  can  be 
changed  in  the  OPEN  statement  in  the  subroutine  transf.f.  This  is  a  large  Unary 
file  and  cannot  be  written  to  the  screen. 
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PARAMETER  FILE  STRUCTURE 

linel*  RED,  which  is  the  reduction  velocity  in  Idkxneters  per  second.  If  RED  is  set 
to  zero,  then  there  is  no  velocity  rediKtion.  A  lumzero  reduction  velocity  gives 
plots  that  have  plotted  times  (P)  equal  to  the  true  travel  time  (T)  less  tiie  (voduct 
of  the  reduction  velocity  (RED)  and  the  offset  (X),  or  P  s:  T  •  (RED  *  X).  A  good 
choice  for  tins  value  is  the  velocity  of  the  fastest  refracting  layer  in  the  model  so 
tiua  the  refractions  from  that  layer  will  iqjpear  flat 

Line  2:  TLAG,  aluch  is  the  tune  lag  in  seconds,  to  be  used  in  the  computed 
synthetics  if  you  want  a  time  offset 

line  3:  JF,  which  is  the  number  of  frequencies  to  be  used  from  the  input  response  file 
("synin**).  If  using  all  the  frequencies  present  in  the  response  file  (output  oi 
PNSEB),  set  JFM). 

Line  4:  NF,  which  denotes  the  number  oi  complex  numbers  used  in  the  fit  NF  must 
be  an  integral  power  of  2  greater  than  or  equal  to  JF.  There  are  2*NF  time 
sanq)les  per  trace  on  the  synthetics.  Normally,  NF  is  set  equal  to  JF,  but  a  larger 
value  can  be  used  to  give  smoother  wave  fimns. 

line  5:  IFLT,  which  is  1  if  earth  flattening  is  used  and  0  otherwise. 

line  6:  IVSP,  which  is  set  to  1  for  VSPs  (NOR  >  1). 

line  7:  INRES,  which  is  set  to  1  for  an  automatic  band  pass  filter.  This  filta  does 
not  preserve  the  anqtlitudes.  The  filters  available  with  lines  13  and  14  are 
preferable. 

line  8:  DPSOURCE,  which  is  set  to  5  for  a  bubble  source  function  taken  from 
Spudich  and  Qtcutt  (1980;  1987)  which  is  not  calibrated  for  amplitude;  set  to  1 
for  Schmidt's  #1  pulse  (Schmidt  1988),  whidi  has  an  amplitude  twice  that  of  a 
sinq)le  iirq)ulse.  Set  to  0  for  no  source,  but  the  calibrated  band-limited  response. 
Note  that  the  calibrated  amplitudes  are  only  for  the  displacement  motion  and  are 
relative  to  the  response  at  1  m  in  any  direction  for  the  pressure  refuse  in  a  liquid, 
1  m  directly  below  a  vertical  iiiq)ulse  source  (M1,M2,M3  =  1,0,0  in  the  "resp" 
file),  and  1  m  horizontally  away  from  a  horizontal  iiiq)ulse  source  (M1,M2,M3  = 
0,0,1).  Note  that  an  impulse  source  of  M1,M2,M3  =  1,0,1  in  an  elastic  medium 
also  gives  a  unit  response  at  1  m  eitiier  directly  horizontally  or  directly  vertically. 

line  9:  IWIND,  which  specifies  the  frequency  window: 

0  s  a  rectangular  window  (strongest  sidelobes) 

1  a  a  triangular  window 

2  s  a  Hanning  window  (weak  sidelobes  and  the  preferred  window) 

3  s  another  Hamiing  window 

4  a  a  Hamming  window 

5  a  a  Nuttall  window  (weakest  sitte  lobes  but  widest  central  peak) 


14 


Line  10:  IFl’ '  AD,  is  set  to  1  tt)  write  the  file  heado’  to  the  file  "pnsyn.out" 
This  gives  diagnostics  for  debugging. 

Line  11:  IFRHEAD,  winch  is  set  to  1  to  write  eadi  xecoid  header  to  die  fiie 
”pnsyn.ouL"  This  gives  diagnostics  for  ddnigging. 

Line  12:  IDVA,  is  set  to  1  to  calculate  di^lacements,  2  for  velocity,  and  3  for 
accekratkxL  Since  diqdacements  have  a  siiqder  pulse,  they  can  ht  used  for 
hydrofdione  pressure  ^thedcs  since  the  uniform  media  (water)  will  not  diange 
tte  tdadve  amplitudes.  The  siiqde  diqilacement  pulse,  is  also  easier  to  convcdve 
with  a  known  source. 

Line  13:  RUTEl,  wfaidi  is  the  low  side  oi  a  mute  window  in  sectuids.  If  a  mute  is 
applied,  dien  all  values  before  RUTEl  are  set  equal  to  die  value  at  time  RUTEl. 

Line  14:  RUTE2,  vriiidi  is  the  end  of  the  mute  window  in  seconds.  If  RUTE2  is  set 
to  0.0,  then  there  is  no  mute  iqiplied.  If  a  mute  is  tqiplied,  then  all  values  after 
RUTE2  are  set  equal  to  the  value  at  time  RUTE2.  This  mute  is  designed  to 
eliminate  arti£EK:ts  dutt  sometimes  appear  befne  and  after  the  arrivals  of  interest 

Line  ISa:  (only  if  IFSOURCE  (line  8)  is  1:)  WEIGHT,0.,0.,0.,0.,0.,0. 

WEIGHT  is  die  cmtral  frequency,  which  should  be  about  one-third  of  the 
maximum  frequency.  The  rest  of  the  variables  are  not  used  for  the  SAFARI 
source.  A  lectangi^  window  (IWIND  =  0)  should  be  used  with  diis  source. 
The  frdlowing  zeroes  are  durmny  values,  since  line  i3a  is  read  by  the  same 
command  as  13b. 

Line  15b:  (only  if  IFSOURCE  (line  8)  is  5:)  WEIGHT,  CONV,  DEPTH,  RANGE, 
BDEC,  BRAT,  DRAT,  where: 

WEIGHT  is  the  weight  of  an  explosive  charge  in  kilograms  or  the  pressure  times 
volume  in  pounds  per  square  inch  and  cubic  inches  of  an  airgun  for  source 
numbers. 

CONV  is  the  ctmversion  factor  for  different  sources 

s  1.0  fOT  TNT  in  kilograms, 

0.639  fOT  TOVEX  in  Idlograms, 

s  3884330.0  for  airguns  in  pounds  per  square  inch  and  cubic  inches. 

DEPTH  is  the  depth  of  the  shot  below  the  water  surface  in  meters. 

RANGE  is  die  range  from  the  source  in  meters  frx’ unit  an^ilitude.  This  number  is 
purely  arbitrary  and  does  not  affect  the  wave  fcmn  or  its  arx^litude. 

BDEC  is  the  bubble  decay  speed.  A  larger  number  makes  the  peaks  sharper,  and 
10.0  w<nks  well. 

BRAT  is  the  bubble  ratio  or  the  amplitude  ratio  between  one  peak  and  the 
preceeding  peak,  and  1.5  to  2.0  works  well. 

DRAT  is  the  dilatation  ratio,  how  much  smaller  the  amplitude  for  the  negative 
pressure  parts  of  the  waves  are,  and  3.0  to  5.0  works  weU.  The  best  loQiws  for 
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BIXEC,  BRAT,  and  IMtAT  wfll  depcad  od  eadi  users  situatian  and  can  be 
delemuned  Bom  experimentation. 

line  16:  IFNOISE,  NFILTN,  where: 

IFNOISE  is  1  to  add  mnse  to  the  seisnx^ram  and  0  to  not  add  noise. 

NFILTN  is  the  nundier  filten  used  on  the  noise  and  can  be  anywhere  from  0  to 
10.  Filtering  die  noise  changes  its  qjectrum  and  can  be  used  to  make  the  noise 
match  the  characteristics  oi  the  receiving  system,  die  environment,  or  whatevo^ 
else  may  be  the  source  of  the  simulated  noise. 

Line  17:  AMFN,  v^ch  is  die  anqiUtude  of  the  noise.  This  is  an  absolute  value,  is  not 
reladve  to  the  signal  strength  wiudi  depoids  on  the  range  and  the  dasdc  model 
used,  and  may  be  changed  by  filtering.  To  get  a  specific  signal  to  nmse  ratio, 
several  trials  may  be  necessary  using  inflection  of  die  seisnograms  to  determine 
the  current  signal  to  noise  rado. 

Lines  18  and  19  are  filter  specification  cards  used  with  the  noise  qption  and  are 
repeated  depending  on  number  of  filters  used  as  determined  by  NFILTN  and 
always  fipear  in  pairs.  If  no  filters  are  to  be  used,  then  these  lines  must  not  be 
present,  not  even  carriage  returns. 

Line  18:  M  for  a  minimum-phase  type  filter  and  Z  for  a  zero-phase  type  filter. 
Miiiimum-phase  filters  are  fdiysically  realizable  and  are  good  for  high-cut  filters. 
Zero-phase  filters  produce  no  i^ase  distcwtion  and  are  good  for  low-cut  filters,  but 
are  not  physically  realizable. 

Line  19:  IFILTYP.IDB,  FI,  F2,  where: 

IFILTYP  is  the  filter  type.  Use  1  for  high-cut  Butterworth,  2  for  low-cut 
Butterwordi,  and  3  for  notch. 

IDB  is  the  slqie  in  dB/octave  and  can  be  any  integer.  A  slope  of  6  dB/octave  or 
less  gives  a  gende  filter,  while  a  slope  of  24  dB/octave  gives  a  hard  filter. 

FI  is  the  comer  fiequency  in  hertz.  For  a  notch  filter,  diis  is  the  lower  comer 
fiequency. 

F2  is  the  upper  ctmier  frequency  of  a  notch  filter  in  hertz.  This  variable  is  not 
used  for  other  filter  types. 

Lines  20  and  21  are  filter  specification  cards  and  are  repeated  depending  on  number  of 
filters  used,  up  to  a  maximum  of  10  and  appear  in  pairs  in  the  identical  fcmnat  as 
lines  18  and  19.  If  no  filters  are  to  be  used,  then  these  lines  must  not  be  present, 
not  even  carriage  returns.  The  program  reads  lines  until  none  are  left  A  zero- 
phase,  low-cut  filter  with  a  ctHner  fiequency  of  10  to  20%  of  the  maximuni 
fiequency  calculated  (NW/TSEC  firom  line  3  of  the  "resp"  file)  is  hi^y 
recommended  to  remove  many  of  the  low-fiequency  p  integration  artifacts.  The 
automatic  filter  on  line  7  will  do  this  but  it  does  not  preserve  the  amplitudes. 
These  artifacts  are  also  gready  reduced  by  increasing  the  number  of  fiequencies 
calculated  in  PNSN. 


OUTPUT  OF  THE  PROGRAM 


There  are  three  output  files  the  program: 

1)  The  output  file,  u^uch  gives  such  conq>utational  details  as  the  beguming  range,  the 

fiiul  range,  range  increment,  cpu  time,  etc.  The  name  of  diis  file  is  "fcHiS." 

2)  The  synthetics  file.  This  is  the  desired  file  which  can  be  plotted  to  generate  synthetic 

seismograms.  The  name  of  diis  file  at  present  defaults  to  "fOTtSO"  and  is  a  large 
ASCQ  file  which  can  be  written  to  the  terminal  screra. 

3)  Further  diagnostic  output  i^pears  in  "pnsyn.ouL” 

NOTE:  All  the  output  files  are  in  ASQL 

STRUCTURE  OF  THE  SYNTHETICS  FILE  rfortaO”) 

Line  1:  NL.  KW,  NP,  NT;  in  (1X,2I5,2I8)  format  where: 

NL  is  the  number  (flayers. 

KW  is  the  number  of  (con^lex)  fiequencies  calculated. 

NP  is  the  number  of  complex  frequencies  used  in  the  FFT.  If  NP  >  KW,  then  the 
FFT  series  is  padded  with  zeroes  before  the  transform  is  made. 

NT  is  the  number  oi  time  samples  in  the  time  series. 

Line  2;  BW,  FW,  ILAG,  DECAY,  BEGP,  FINP,  PWl,  PW4,  PW2,  PW3;  in 
(1X,10E14.7)  format  where: 

BW  is  die  starting  frequency  in  hertz. 

FW  is  the  ending  frequency  in  hertz. 

TLAG  is  the  time  lag  in  seconds. 

DECAY  is  the  decay  factor  over  the  time  series  length. 

BEGP  is  the  beginning  ray  parameter. 

FINP  is  the  final  ray  parameter. 

PWl  is  the  lower  edge  of  the  ray  parameter  (p)  window. 

PW4  is  the  upper  edge  of  die  p  window. 

PW2  is  the  start  of  the  taper  for  the  lower  part  of  the  p  window. 

PW3  is  the  start  of  the  taper  for  the  upper  part  of  the  p  window. 

Line  3:  T(I),  CL(D,  CTO).  QL(D,  QT(I).  RHO(I);  in  (1X,6E14.7)  format  and 
appears  once  for  each  layer. 

T(I)  is  die  thickness  of  layo*!  in  kilometers. 

CL(I)  is  the  compressional  wave  velocity  for  layer  I  in  kilometers  per  second. 

CT(I)  is  the  shear  wave  velocity  for  layo*  I  in  kilometers  per  second. 

QL(I)  is  the  ccmqnessitmal  Q  for  layer  L 
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QIXI)  is  die  shear  Q  for  layer  L 

RHO(I).is  die  (hmsity  for  layer  I  in  gnons  per  cubic  centimeter. 

Line  4:  NFILT,  NONPHASE.  VSP,  LS(XJR.  NOR*  NX.  ZS.  RV.  BX.  FX.  DX  in 
(1XJ2.2L1.3I3.5E14.7)  fonnat  where: 

NFILT  is  die  number  offilters(nnxinaim(tf  10). 

MINPHASE  is  tine  if  the  filters  are  minimum  phase,  else  false. 

VSP  is  true  if  there  is  more  than  <we  receiver  dqMfa. 

LSOUR  is  the  layer  number  that  contains  the  source. 

NOR  is  the  number  of  receiver  dqiths.  (1  for  normal  seismograms) 

NX  is  the  number  of  receiver  offiKts. 

ZS  is  the  depth  of  die  receiver  in  kilometers. 

RV  is  the  reducing  velocity  in  kilometers  per  second. 

BX  is  the  offset  of  the  first  receiver. 

FX  is  the  offset  of  the  last  receiver. 

DX  is  the  difference  in  offset  between  receivers. 

lines*.  IFILTYP(10).IDB(10)  in  (1X.10I1.10I3)  fonnat  where: 

IFILTYP  are  the  filter  types  for  all  10  filters;  1  for  hi|^  cut.  2  for  low  cut.  3  for 
notch. 

IDB  are  die  filter  slqies  in  dB/octave  for  all  10  filten. 

Lined:  Fl(lO)  in  (1X,10F10.4)  fonnat  where: 

FI  are  the  frequencies  in  hertz  for  all  10  filters.  For  a  notch  filter,  this  is  the  lower 
finquency. 

Line  7:  F2(10)  in  (1X,1(F10.4)  format  where: 

F2  are  the  upper  frequency  in  hertz  for  notch  filters. 

Lines  8  and  9  are  repeated  in  sequence  cmce  for  each  trace  (or  offset)  in  the  computed 
synthetics. 

Lines;  ISS,  NT.  DT.  X.  Z;  in  (1X.218,3E14.7)  FORMAT  where: 

ISS  is  a  flag  for  the  end  of  the  file.  It  is  0  for  every  trace  excqit  for  the  last  trace 
when  it  is  1. 

NT  is  the  nundier  of  time  samples. 

DT  is  die  sanqiling  interval  in  seconds. 

X  is  the  range  in  kilometers  for  this  trace. 

Z  is  die  dqpdi  in  kilometers  for  diis  recdver. 
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line  9  is  rq>eated  NT/8  times  fw  each  trace  and  contains  the  time  series  values,  8  at  a 
time,  oi  the  trace  in  8(1X,E14.7)  fonnat  The  smplimde  of  the  time  ^esies 
relative  to  a  direct  arrival  at  one  meter. 

RUNNING  THE  VERSION  SYNSPEC 

SYNSFEC  uses  the  sanlh  format  fw  the  "param"  file.  In  £u»,  all  the  subroutines 
are  die  same  except  for  piisyn.f,  which  bectxries  synspec.f.  The  makefile  and  the  runfile 
for  SYNSPEC  are  the  same  as  for  PNSYN  except  that  "syn^iec"  is  used  in  place  oi 
"pns3fn." 


I  ,  .  lu  Mill,  1..  ,1.111 . III.  I  1.11.  ,  L.J.  .  1. 


SAMPLE  FILES: 

Ftdlowing  are,  in  older,  a  "inakefile"^  PNS  YN,  a  shell  script  file  to  nin  PNS  YN,  and 
two  "paiam**  files  id  show  cases  with  and  widuNit  sources  and  films. 

#  Makefile  for  pnsyn 

#  Makefiles  are  used  In  UNIX  based  systems  and  are 

#  system  dependent.  The  best  way  for  you  to  make  such  a 

#  file  Is  with  the  "makemake”  command  and  then  edit  It  to 

#  Include  the  options  that  you  want.  This  makefile  Is  for 

#  a  CONVEX  system  and  Includes  some  vectorlzatlon  and 

#  optimization  options. 

OBJ  =\ 

pnsyn. o  parmsln.o  transf.o  sources. o  pulsel.o  filter. o  bufln.o  \ 
fft.o  get secs. o 

FC  =  fc  -vn 

OPT  =  -01  -cfc 

PROF  = 

LINKPROF  = 

LIB  =  -Ivecllb 

.  f  .o: 

$(FC)  $(OPT)  $(PROF)  -c  $★.! 
pnsyn. 01:  $ (OBJ) 

$(FC)  $ (LINKPROF)  -01  -o  pnsyn. 01  $ (OBJ)  $ (LIB) 

get secs. o: 

$ (FC)  $ (PROF)  -c  getsecs.f 
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#  Makefile  for  synspec 

OBJ  -\ 

pnsyn.o  parmsln.o  transf.o  sources. o  filter. o  bufln.o  bufout.oX 
fft.o  get secs. o 

FC  -  fc  -vn 

OPT  »  -01  -cfc 

PROF  » 

LINKPROF  » 

LIB  =*  -Ivecllb 

.  f  .o: 

$(FC)  $(0PT)  $(PR0F)  -c  $^.f 

synspec. 01:  $ (OBJ) 

$(FC)  $ (LINKPROF)  -01  -o  synspec. 01  $ (OBJ)  $ (LIB) 

getsecs . o : 

$(FC)  $(PR0F)  -c  getsecs. f 
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The  next  file  is  t  diell  soipt  (UNIX)  for  running  the  program  ^SYN.  It  assigns 
logical  nmobers  to  files  and  directs  error  li^g  to  "imsn.oiu."  Note  diat  the  data  file 
produced  by  pnsn  (deiqi.  denh,  or  denv)  has  been  renaii^  "synin." 

#  11  is  a  parameter  Input  file 

setenv  forOll  /lo2way/llndwall/pn/syn/pararo 

#  20  Is  the  w-x  domain  synthetics  from  pnsn 
setenv  FOR020  /lo2way/llndwall/pn/synln 
unallas  rm 

rm  syn.out 

/lo2way/llndwall/pn/syn/pnsyn.01  >S!  syn.out 

#  fort. 30  Is  the  time  domain  output  file 

#  fort. 3  Is  the  diagnostic  output  file 
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The  next  two  files  are  sample  "param"  files.  The  first  is  for  PNSYN  and  includes  a 
filter  and  a  source  function  line.  The  source  function  is  for  the  bubble  pulse  fixnn  27  kg 
of  TOVEX.  This  file  is  set  for  a  PNSYN  input  file  ("synin")  calculated  with  512 
fiequencies.  NF  is  set  to  1024  rather  than  512  so  that  the  waveforms  will  appear 
smoother  in  the  plots.  A  2-Hz,  high-pass  filter  is  used  to  remove  some  of  the  voy  low 
frequency  artifacts.  The  second  "param"  file  is  for  a  SYNSPEC  input  file.  This  file  could 
also  include  filters  and  a  source  but  does  not 


6.0 

0. 

512 

1024 

0 

0 

0 

5 

2 

0 

0 

1 


1.1 

0.0 

27.3.0,0.639,78 


!  reduction  velocity  in  km/s 
!  time  lag  in  sec. 

!  number  of  frequencies  to  use 
!  number  of  points  to  use  in  the  fft 
!  set  to  1  if  earth  flatening  is  used 
!  set  to  1  for  VSPs 

!  set  to  1  for  automatic  band  pass  filter 
!  specifies  source,  0  is  no  source 
!  specifies  frequency  window,  2  is  banning 
!  set  to  1  to  write  file  header 
!  set  to  1  to  write  record  headers 
!  1-displacement,  2-velocity,  3-acceleration 
!  mute  window  1,  seconds 
!  mute  window  2,  if  0.0,  then  no  mute 
.,10000. ,10. ,1.5, 5. 


Z 


2,12,2. 


0.0 

0. 

512 

512 

0 

0 

0 

0 

0 

0 

0 

1 

1.1 


0.0 


reduction  velocity  (km/s) 
time  lag 

length  of  time  series  to  read 
length  of  fft  time  series 
1  to  use  earth  flatening 
1  for  VSPs 

1  for  band  pass  filter 

source  type,  0  for  no  source 

0  for  amplitude,  1  for  power  spectra 

set  to  1  to  write  file  header 

set  to  1  to  write  record  headers 

1-displacement,  2-velocity,  3-acceleration 

mute  window  1,  seconds 

mute  window  2,  if  0.0,  then  no  mute 
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VSPLT 


The  fiks  from  PNSYN  or  SYNSPEC  are  ASCII  files  and  can  be  transported  to  many 
other  computers.  Once  the  files  are  on  the  cmiqjuter  to  do  the  plotting,  tte  data  nuist  be 
converted  to  a  binaiy  file  with  the  program  CONVRT.  CONVRT  requires  no  user 
parameter  file.  It  te^  the  "fortSO"  file  and  a  title  file  "synhead"  and  writes  a  binary 
"fortSr  data  file  and  a  binaiy  "synhead"  fite,  vMch  contains  smne  heatter  infonnation. 
After  running  CONVRT,  the  plotting  is  done  by  running  VSPLT.  Several  lines  in  tiie 
FORTRAN  program  VSPLT  are  plotter  dependent  and  plotting  software  dependent  and 
may  need  to  be  changed  for  the  users  systeoL  This  code  has  been  run  on  many  madiines 
and  the  iiutializing  calls  for  several  types  of  plotters  are  still  in  the  code  as  comment  lines. 

The  package  VSPLT  cmisists  of  the  main  program  and  four  subroutines: 

1)  vsplLf,  the  main  module. 

2)  The  subroutine  axis.f,  which  plots  the  axis  fOT  the  section. 

3)  The  subroutine  plottr.f,  which  plots  the  time  series  traces. 

4)  The  subroutine  step.f. 

5)  The  sulnoutine  modplLf,  which  plots  the  geophysical  model  before  the  time  series. 

On  UNIX  systems,  the  omipiladon  and  linking  is  most  ctmveniently  done  with  a 
makefile.  A  sample  "makefile”  is  included  with  the  set  of  sample  iiq>ut  files. 

ROSECON  reads  a  file  from  PNSYN  and  writes  an  ASCII  file  in  ROSE  fenrnat 
(LaTtaille  and  Dorman  1982;  LaTraille  et  al.  1983)  which  can  be  plotted  by  several 
programs  desrgned  for  that  purpose.  While  CONVRT  has  no  input  file,  ROSECON  does 
to  specify  certain  variable  parameters. 

INPUTS  TO  VSPLT 

The  inputs  to  the  program  VSPLT  are  in  one  file,  "paramS".  The  name  can  be 
changed  in  the  shell  script  "runplot"  In  this  input  file,  the  names  of  integer  variables  begin 
with  the  letters  I-N.  All  others  variables  are  real  numbers  whidi  must  have  a  decimal 
pmnL 

The  "paiamS"  file  contains  the  following  in  sequence: 

Line  1:  IPLOTYP  specifies  the  receiver  type;  1  is  for  pressure,  2  is  for  horizontal  and 
3  is  for  vertical  This  is  used  for  information  only  and  does  not  affect  the  plot 

line  2:  RMIN  is  the  miniimim  range  on  the  plot  or  the  ori^n  of  the  x  axis  in 
kilometers. 

Line  3:  RAMX  is  the  maximum  range  on  the  plot  in  kilometers.  Occasionally  the 
program  will  plot  to  a  slightly  greater  range. 
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line  4:  RSF  is  the  range  scaling  factor.  This  is  in  units  of  inches  per  100  km. 

Line  5:  TMIN  is  the  miniinum  time  on  the  plot  in  seconds.  This  can  either  be  the 
origin  or  the  top  of  the  plot,  depending  on  the  sign  of  TSF. 

Line  6:  NET  is  the  number  of  seconds  to  plot,  an  integer.  The  maximum  time  on  the 
plotisTMIN  +  NET. 

Line  7:  TSF  is  the  time  scaling  factiOT  with  units  of  inches  per  second.  When  diis 
factor  has  a  negative  sign,  TMIN  will  be  at  the  tc^  of  the  plot  and  time  will 
increase  downward. 

Line  8:  AMP  is  the  gain  vaiiaUe.  All  data  values  are  multiplied  by  diis  factor  before 
plotted. 

Line  9:  DIR  is  the  depth  to  the  first  reflector  used  by  the  mute.  When  this  value  is  0, 
no  mute  is  used. 

Line  10:  VIR  is  the  velocity  overlying  the  first  reflector  and  is  used  to  mute  the  noise 
in  front  of  the  first  reflector  hyperbola. 

Line  11:  GPS  is  the  gain  per  second  after  the  first  reflector.  When  GPS  =  0.0,  no 
dme-vaiying  gain  is  used 

Line  12:  ALPHA  is  the  exponent  used  for  a  range  dependent  gain.  Whm  ALPHA  = 
0.0,  no  range  dependent  gain  is  used.  Also,  be  sure  to  set  ALPHA  to  0.0  if  there  is 
a  record  at  zero  offset. 

Line  13:  XXO  is  the  reference  distance  for  range  dependent  gain.  The  default  value  is 
1  km. 

Line  14:  ISCALE  selects  a  trace  dependent  amplitude  scaling.  If  ISCALE  =  1,  then 
the  maximum  amplitude  of  each  trace  is  scaled  to  the  maximum  amplitude  of  the 
first  trace.  If  ISCALE  =  2,  then  each  trace  is  scaled  to  its  own  maximum 
aiiq)litude.  If  ISCALE  =  3,  then  no  trace  scaling  is  done  and  each  trace  is 
multiplied  only  by  the  AMP  factor  and  tiie  range  scaling  factor. 

Line  IS:  lOPT  is  the  switch  for  shading.  Setting  lOPT  =  1  gives  a  variable  area  plot 
(witii  shading)  while  setting  lOFT  =  0  gives  a  wiggle  line  plot 

Line  16:  ISGN  selects  which  side  is  shaded  when  using  the  variable  area  plot  option 
(lOPT  s  1).  1  shades  the  +  side  while  *1  shades  the  -  side. 

Line  17:  MODPLOT  selects  the  option  of  plotting  the  gecq)hysical  model  in  firont  of 
the  record  section.  MODPLOT  =  1  plots  the  nxxlel  while  MODPLOT  =  0  does 
not 

Line  18:  ICLIP  selects  whether  or  not  to  clip  the  traces.  This  is  highly  recommended. 
Set  ICUP  s  1  to  clip  the  traces. 

Line  19:  CLAMP  which  is  the  maximum  amplitude  from  positive  peak  to  negative 
peak  in  inches  allowed  by  the  clipper. 

Line  20:  TDELAY  adds  TDELA Y  rime  in  seconds  to  the  front  of  the  traces. 
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Line  21:  JDB  selects  the  logarithmic  plotting  q>tion.  JDB  »  1  dien  the  (riot  is  in 
decibels,  if  JDB  »  0,  then  the  traces  are  i^ott^  lineariy. 

Line  22:  DBOFF  is  the  DC  offset  added  to  each  trace  when  plotted  in  decibels. 

Line  23:  ISTART  is  the  first  trace  to  be  plotted. 

line  24:  UUMP  is  the  interval  between  traces  to  plot.  Every  UUMPth  trace  is 
plotted. 

Line  25:  NTOPLOT  is  the  number  of  tracres  to  plot  NTOPLOT  s  0  plots  all  traces. 

OUTPUT  OF  THE  PROGRAM 

The  program  produces  a  vector  file  which  must  be  converted  to  something  readable  by 
your  madiine.  This  may  be  a  system  software  command  like  "raster,"  which  makes  a 
raster  rile  for  a  dot-matrix  plotter. 

RUNNING  ROSECON 

ROSECON  reads  the  same  data  file  as  CONVRT  but  requires  a  parameter  file 
"roseconin."  The  first  six  lines  of  the  "roseconin"  file  have  one  parameter  each,  which  are, 
in  mden  instrument  number,  event  number,  channel  number,  instrument  ID  number, 
source  type,  and  receiver  type.  The  last  two  lines  are  shot  time  and  first  sample  time  with 
sue  parameters  on  each  line  being,  in  order  year,  month,  day,  hour,  minute,  second.  A 
"roseconin"  file  is  included  with  the  sanq)le  files. 

OTHER  PLOTTING  PROGRAMS 

Other  plotting  programs  will  inevitably  be  written  as  printers  and  software  are  updated 
and  may  not  use  CONVRT.  Contact  other  users  to  see  what  is  currently  available. 

SAMPLE  FILES: 

On  the  following  pages  are,  in  order,  a  "roal^e"  for  VSPLT,  a  shell  script  file  to  nm 
VSPLT,  a  "paiamS"  file,  makefiles  for  CONVRT  and  ROSECON,  shell  script  files  to  run 
CONVRT  and  ROSECON,  and  a  "roseconin"  rile. 
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« 


Makefile  for  vsplt 


OBJ  »\ 

vsplt. o  axls.o  plottr.o  step.o  modplt.o 
FC  »  fc  “vn 
OPT  =  -01 
PROF  = 

LINKPROF  - 

LIB  =  -Ivpltlib  -Iveclib 
.  f  .o: 

$(FC)  $(0PT)  $(PR0F)  -C  $it.f 
vsplt. 01:  $(0BJ) 

$(FC)  $ (LINKPROF)  -01  -O  vsplt. 01  $ (OBJ)  $ (LIB) 


#  Sample  script  to  run  plotting  code.  These  commands  are 

#  mashine  and  software  dependent.  Consult  your  local 

#  computer  gurus 
unalias  rm 

rm  /io2way/lindwall/pn/V9INT.DAT 
rm  /io2way/lindwall/pn/VERSATEC 

#  assign  lun  11  to  plot  parameter  file  (you  generate) 
setenv  FOROll  /io2way/lindwall/pn/plot/param5 

#  assign  model  parameter  file  (lun  11  generated  by  convrt)  to 

#  lun  12 

setenv  FOR012  /io2way/llndwall/pn/synhead 

#  assign  lun  10  to  the  fort.xx  (xx  «  81,82,..)  generated  by 

#  the  convrt  program 

setenv  FOROlO  /io2way/lindwall/pn/fort . 81 
/io2way/lindwall/pn/plot/vsplt .01  >&!  vsplt. out 

« 

#  the  next  two  commands  plot  the  record  section 
raster 

dd  if=VERSATEC  of  *  /dev/pbS  bs=528 
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OU»OMOOO»-*0lN>O 


The  neia  fite  is  t  saovte  "pnaaS”  fik  whidi  qwdfies  all  of  the  pkx  aiwcificatioiis. 

1  !  Iplotyp,  l«pres,  2«horlz,  3«vert 

0.0  !  rmln,  minimum  range  in  km. 

20.0  !  rmax,  maximum  range  in  km. 

20.0  !  rsf,  lnches/100  km 

.0  !  tmln,  minimum  time  on  plot 

!  net,  number  of  seconds  to  plot 
.0  !  tsf,  Inches/sec;  when  <0,  time  Increases  downward 

.0  !  amp 

.0  !  dir,  depth  to  first  reflector 

.  !  vlr,  velocity  overlying  first  reflector 

.  !  gps,  gain  per  second  after  first  reflector 

.0  !  alpha,  range  exponent,  use  0  if  there  is  an  x»0  trace 

.  !  xxO 

!  Iscale,  l»flrst  trace,  2«max  of  each,  3«not  scaled 
!  lopt,  set  to  one  (1)  for  shading  (variable  area) 

-1  !  Isgn,  1  shades  the  side,  -1  shades  the  -  side. 

0  !  modplt,  set  to  1  to  get  plot  of  model. 

1  !  Icllp,  if  »  1,  clips  at  "clamp"  Inches 

2.0  !  clamp 

0.0  !  tdelay,  delays  record  by  tdelay  seconds 

0  !  jdb,  if  1,  plot  in  dB  Instead  of  linear  amplitude 

-00.  !  dboff,  dc  offset  for  dB  values 

1  !  istart,  first  trace  to  plot 

1  !  ijump,  plot  every  "ijump"th  trace 

0  !  ntoplot,  number  of  traces  to  plot.  0  plots  all 
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« 


Makefile  for  convrt 


OBJ  »\ 
convrt . o 

FC  *  fc  -vn 

OPT  =  -01 

PROF  = 

LINKPROF  - 

LIB  *  -Iveclib 

.  f  .o: 

$(FC)  $(OPT)  $(PROF)  -c  $^.f 

convrt. 01:  $ (OBJ) 

$(FC)  $ (LINKPROF)  -01  -o  convrt. 01  $ (OBJ)  $ (LIB) 

#  Makefile  for  rosecon 

OBJ  =\ 
convrt . o 

FC  =  fc  -vn 

OPT  =  -01 

PROF  = 

LINKPROF  » 

LIB  «  -Iveclib 

.f  .o: 

$(FC)  $(OPT)  $(PROF)  -c  $^.f 

rosecon . 01 :  $ (OBJ) 

$(FC)  $ (LINKPROF)  -01  -o  rosecon.Ol  $ (OBJ)  $ (LIB) 
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The  next  two  files  are  sample  shell  scripts  for  running  CX)NVRT  and  ROSEOCW. 


» 

unallas  rm 

rm  convrt.out  synhead  fort. 81  fort. 82 

« 

#  10  Is  the  Input  synthetics 
setenv  FOROlO  fort. 30 

#  11  Is  the  heading  output  file 
setenv  FOROll  synhead 

#  12  Is  the  title  Input  file 
setenv  FOR012  stltle 

#  following  are  files  of  traces,  one  file  per  offset  for  VSP 

#  or  one  file  per  receiver  depth  otherwise.  Start  with  unit 

#  81,  but  don't  exceed  unit  98. 

# 

I  setenv  FOR081  synrec 

« 

/lo2way/llndwall/pn/con/convrt .01  >4i  convrt.out 


« 

unallas  rm 

rm  convrt.out  synhead  fort. 81  fort. 82 

« 

#  10  Is  the  Input  synthetics 
setenv  FOROlO  fort. 30 

#  11  Is  the  heading  output  file 
setenv  FOROll  synhead 

#  12  Is  the  title  Input  file 
setenv  FOR012  stltle 

#  21  Is  an  Input  file  with  rose  format  Information 
setenv  FOR021  /lo2way/llndwall/pn/con/roseconln 

#  following  are  flies  of  traces,  one  file  per  offset  for  VSP 

#  or  one  file  per  receiver  depth  otherwise.  Start  with  unit 

#  81,  but  don't  exceed  unit  98. 

« 

#  setenv  FOR081  synrec 

» 

/lo2way/llndwall/pn/con/rosecon.01  >&!  rose. out 
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This  is  a  "loseoonin''  file. 


1  !  Instrvunent  number 

1  !  event  number 

1  !  channel  number 

1  !  Instrument  ID  number 

2  !  1  for  earthquake,  2  for  shot 

4  !  1  vertical;  2  radial;  4  hydrophone 

90,1,1,0,0,0  !y,m,d,h,m, s  for  shot  and  first  sample  time 
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EXAMPLES 


FoUowiiig  are  sevtxal  examples  calculated  and  plotted  to  show  what  die  synthetic 
seismognuns  should  look  like.  Each  exan^le  consists  of  a  "model"  file  m  the  upper  left 
corner,  a  "re^"  file  in  the  iqiper  right  comer,  followed  by  a  "param"  file  ftx-  nmning 
PNSYN,  and  dien  the  seismic  record  section.  The  examples  shown  are: 

1)  A  horizontal,  near  aea>sorftioe  line  w^  a  simple  subbottom.  The  strong  airivals  at 
about  3  to  5  km  (^aet  and  9  s  reduced  travel  time  are  artifacts  the  calculation. 

2)  A  horizontal,  near  sea-surfKe  line  with  a  simplified  ocean  crustal  model 

3)  A  horizontal,  near  sea-surface  line  with  m  ocean  crustal  model  induding  gradients. 
The  record  section  is  plotted  with  a  different  rethiction  velocity  than  those  above  to  better 
show  die  long  ran^  reftaction  arrivals  that  are  typical  in  marine  seismic  refraction  lines. 

4)  A  line  of  near  sea-surface  shots  to  an  ocean  bottom  seismometer  (OBS).  This 
inclu^  a  record  section  eadi  for  die  pressure,  horizontal,  and  vertical  conqionents.  The 
feature  seen  at  20  km  and  about  2  s  is  a  wnyiaround  artifact 

5)  A  vertical  jaeignic  profile  (VSP).  The  receivers  are  jdaced  at  the  top  of  layers  3 
through  12  (the  variaUe  LOBS  in  the  "resp"  fileX  and  the  source  is  at  the  top  of  layer  13 
at  a  range  of  SOOm  (variable  BX).  The  program  wiU  not  work  conecdy  if  the  source  is  in 
die  same  layer  as  one  of  the  tecdvets. 

The  surface  lines  (exarrqiles  1,  2,  and  3)  can  be  either  a  single  shot  to  a  line  of 
receivers  or  a  line  of  shots  to  a  single  receiver.  In  these  calculations,  the  problems  are 
reciprocal  likewise,  the  OBS  line  could  alternatively  be  thought  of  as  a  sin^  surface 
shot  to  a  line  a(  OBSs. 
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0.07, 50. ,0,3, 0,0 

1. 5. 5000. . 0..0.. 1.03. 0.01 
1. 5,  5000. ,  0.,  0.,  1.03,  0.005 

1.5.5000.. 0..0..1.03.3.985 

1.8. 100. . 0.6. 6. .1.4. 0.5 

4. 5. 100. . 2. 2,  50. ,2. 3,  50.0 


: MODEL  HSP  1 
:PVH 

256,1,10.0,0.8 

2 

1.0,20.0,1.0 
0., 0.1, 0.7, 0.8 
0.005 
D 

0.,0.,0. 

1., 0.0,1. 


!nw, nor, tsec,  p2w 
Hobs  (nor) 

!bx,  fx,dx 

! pwl , pw2 , pw3 , pw4 

Idpmax 

! otype 

!h(l) ,h(2) 

!  ml ,  m2 ,  m3 


4.0 

0.0 

256 

512 

0 

0 

0 

0 

2 

0 

0 

1 


1.1 

0.0 

0,0 

0.020 

Z 

2,12,2. 

M 

1,12,13. 


!  reduction  velocity  in  km/s 
!  time  leg  in  sec. 

!  number  of  frequencies  to  use 

!  number  of  points  to  use  in  the  fft  (power  of  2) 

!  set  to  1  if  earth  flatening  is  used 
!  set  to  1  for  VSPs 

!  set  to  1  for  automatic  band  pass  filter 
!  specifies  source,  0  is  no  source 

!  specifies  frequency  window,  0  is  rectangle,  2  is  banning 
!  set  to  1  to  write  file  header 
!  set  to  1  to  write  record  headers 

!  I'displacement,  2-velocity,  3-acceleration  (pressure) 

!  mute  window,  seconds 

!  end  of  mute  window,  if  0,  then  no  mute 
.'  1  for  noise,  number  of  noise  filters 
!  noise  amplitude 

horizontal  profile  #  1 
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0.0'}, so. ,0,3, 0,0 
X.S,SO00.,0.,0., 1.03, 0.01 

1. 5. 5000.. 0..0.. 1.03. 0.005 

1.5.5000.. 0..0..1.03.3.985 

1.8. 100. . 0.6. 6. .1.4. 0.5 

4. 5. 100. . 2. 2. 50. .2. 3. 0.5 

5. 7. 100. . 2. 2. 50. .2. 3. 1.0 

6. 8. 100. . 3. 6. 50. .3. 0.4. 50 
8. 0,200., 4. 5, 100., 4. 5, 50. 00 


:HI»EL  HS?  2 
t  IPVH 

2^6,1,10.0,0.8 

2 

1.0,20.0,1.0 
0., 0.1, 0.7, 0.8 
0.005 
D 

0.,0.,0. 

1., 0.0,1. 


!nw,no£,t9«c,p2w 
!  lobs (nor) 

!bx,  fx,  dx 

!  pwl , pw2 , pw3, pw4 

.'dpfflsx 

!etyp« 

!h(l)  ,h(2) 

!  ml ,  m2 ,  m3 


4.0 

0.0 

256 

512 

0 

0 

0 

0 

2 

0 

0 

1 

1.1 

0.0 

0,0 

0.020 

Z 

2,12,2. 

M 

1,12,13. 


!  reduction  velocity  in  km/s 
!  time  leg  in  eec. 

!  number  of  frequencies  to  use 

!  number  of  points  to  use  in  the  fft  (power  of  2) 

!  set  to  1  if  earth  flatening  is  used 
!  set  to  1  for  VSPs 

!  set  to  1  for  automatic  band  pass  filter 
!  specifies  source,  0  is  no  source 

!  specifies  frequency  window,  0  is  rectangle,  2  is  banning 
!  set  to  1  to  write  file  header 
!  set  to  1  to  write  record  headers 

!  1-displaceraent,  2-veloclty,  3-acceleration  (pressure) 

!  mute  window,  seconds 

!  end  of  mute  window,  if  0,  then  no  mute 
!  1  for  noise,  number  of  noise  filters 
!  noise  anplitude 


horizontal  profile  #  2 


O 

s 


10  15 

offset  in  km 
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0.07, 50. ,0,  3, 0,0 

1. 5. 5000. . 0..0.. 1.03. 0.01 

1. 5. 5000. . 0..0.. 1.03. 0.005 

1.5.5000.. 0..0..1.03.3.985 

1.7. 100. . 0.5. 6. .1.4. 0.1 
0.,3.,0.3 

1.9. 100. . 0.8. 6. .1.4. 0.1 
4. 0,100. ,2. 0,50. ,2. 3, 0.1 
0.,5.,1.2 

5. 7. 100. . 3. 2. 50. .2. 3. 0.2 
0.,5.,2.0 

6. 8. 100. . 3. 6. 50. ,  3. 0,  0. 20 
0.,5.,2.1 

7. 2. 100. . 3. 8. 50. .3. 0.0. 20 
8. 0,200. ,4. 5, 100. ,4. 5, 50. 00 


: MODEL  HSP  2 -GRAD 
:PVH 

256,1,10.0,0.8 

2 

1.0,40.0,1.0 
0., 0.1, 0.7, 0.8 
0.005 
0 

0.,0.,0. 

1., 0.0,1. 


! nw,  nor,  tsec, p2w 

! lobs (nor) 

! bx,  fx,  cbc 

!pul,pw2,pw3,pw4 

idpmax 

!etyp« 

!h(l)  ,h(2) 

!  ml ,  m2 ,  in3 


6.0 

0.0 

256 

512 

0 

0 

0 

0 

2 

0 

0 

1 

1.1 

0.0 

0,0 

0.020 

Z 

2,12,2. 

M 

1,12,13. 


!  reduction  velocity  in  (cin/s 
!  time  lag  in  sec. 

!  number  of  frequencies  to  use 

!  number  of  points  to  use  in  the  fft  (power  of  2) 

!  set  to  1  if  earth  flatening  is  used 
!  set  to  1  for  VSPs 

!  set  to  1  for  automatic  band  pass  filter 
!  specifies  source,  0  is  no  source 

!  specifies  frequency  window,  0  is  rectangle,  2  is  banning 
!  set  to  1  to  write  file  header 
!  set  to  1  to  write  record  headers 

!  l-displacement,  2-velocity,  3-acceleration  (pressure) 

!  mute  window,  seconds 

!  end  of  mute  window,  if  0,  then  no  mute 
!  1  for  noise,  number  of  noise  filters 
!  noise  anq^litude 

HSP  #  2  with  gra(jients 


offset  in  km 
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0.07, 50. ,0,2, 0,0 

1. 5,  5000. ,  0.,0., 1.03,  0.01 

1.5.5000.. 0..0..1.03.3.99 

1.8. 100. . 0.6. 6. .1.4. 0.5 

4. 5. 100. . 2. 2. 50. .2. 3,  0.5 

5. 7. 100. . 2. 2. 50. .2. 3. 1.0 

6. 8. 100. . 3. 6.50. .3. 0.4. 50 
8. 0,200. ,4. 5, 100. ,4. 5, 50. 00 


:OBS  MODEL  2 
:PVH 

256,1,10.0,0.8 

3 

1.0,20.0,1.0 
0., 0.1, 0.7, 0.8 
0.005 
D 

0.,0.,0. 

1., 0.0,1. 


! nw, nor, tsec, p2w 
! lobs (nor) 

! bx,  fx^dx 

!pwl,pw2,pw3,pw4 

! dpmax 

letype 

!h{l)  ,h(2) 

!  ml ,  m2 ,  m3 


6.0 

0.0 

256 

512 

0 

0 

0 

0 

2 

0 

0 

1 

1.1 

0.0 

0,0 

0.020 

z 

2,12,2. 

M 

1,12,13. 


!  reduction  velocity  in  )un/s 
!  time  lag  in  sec. 

!  numiDer  of  frequencies  to  use 

!  number  of  points  to  use  in  the  fft  (power  of  2) 

!  set  to  1  if  earth  flatening  1?  used 
!  set  to  1  for  VSPs 

!  set  to  1  for  automatic  band  pass  filter 
!  specifies  source,  0  is  no  source 

!  specifies  frequency  window,  0  is  rectangle,  2  is  banning 
!  set  to  1  to  write  file  header 
!  set  to  1  to  write  record  headers 

!  1-displacement,  2-velocity,  3-acceleration  (pressure) 

!  mute  window,  seconds 

!  end  of  mute  window,  if  0,  then  no  mute 
!  1  for  noise,  number  of  noise  filters 
!  noise  amplitude 

OBS  -  pressure 


offset  in  km 
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« 


OBS  -  vertical 


5  10 


offset  in  km 


OBS  -  horizontal 


offset  in  km 


39 


0.07,50. , 0, 13, 0, 0 

1. 5.5000. . 0..0.. 1.03. 0.01 

1. 5.5000. . 0..0.. 1.03. 0.99 
1 . 8,-100. ,  0. 6,  €.,1.4,0.05 

1.8. 100..  0.6, €.,1.4, 0.05 

1.8. 100. . 0.6. 6. .1.4. 0.05 

1.8, 100. ,0.6, 6. ,1.4, 0.05 

1.8, 100. ,0.6, 6. ,1.4, 0.05 

1.8, 100. ,0.6, 6. ,1.4, 0.05 

1.8, 100. ,0.6, 6. ,1.4, 0.05 

1.8, 100. ,0.6, 6. ,1.4, 0.05 

1.8, 100. ,0.6, 6. ,1.4, 0.05 

1.8. 100. . 0.6. 6. .1.4. 0.05 

4. 5. 100. . 2. 2,  50. ,2. 3, 0.5 

5. 7. 100. . 2. 2. 50. .2. 3. 1.0 

6. 8. 100. . 3. 6,  50. ,3. 0,4. 50 
8. 0,200. ,4. 5, 100. ,4. 5, 50. 00 


:HODeL  2 
•  pVH 

256,10,10.0,0.8 

3,4,5,6,7,8,9,10,11,12 

0.5, 0.5, 0.5 

0., 0.1, 0.7, 0.8 

0.005 

D 

0.,0.,0. 

1., 0.0,1. 


!nw,nor,tsec,p2w 

! lobs (nor) 

!bx, fx,dx 
!pwl,pw2,pw3,pw4 
!  dpmax 
!«typ« 

!h(l) ,h(2) 

!  ml ,  in2 ,  m3 


0.0 

0.0 

256 

512 

0 

1 

0 

0 

2 

0 

0 

1 


1.1 

0.0 

0,0 

0.020 

Z 

2,12,2. 

M 

1,12,13. 


!  roduction  velocity  in  km/s 
!  time  lag  in  sec. 

!  number  of  frequencies  to  use 

!  number  of  points  to  use  in  the  fft  (power  of  2) 

!  set  to  1  if  earth  flatening  is  used 
!  set  to  1  for  VSPs 

!  set  to  1  for  autOMtic  band  pass  filter 
!  specifies  source,  0  is  no  source 

!  specifies  frequency  window,  0  is  rectangle,  2  is  banning 
!  set  to  1  to  write  file  header 
!  set  to  1  to  write  record  headers 

!  1-displacement,  2-velocity,  3-acceleration  (pressure) 

!  mute  window,  seconds 

!  end  of  mute  window,  if  0,  then  no  mute 
!  1  for  noise,  numlser  of  noise  filters 
!  noise  anplitude 

VSP  -  vertical 
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